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Abstract 

Both IrV and RhV crystalhze in the alrV structure, with a transition to the 
higher symmetry LIq structure at high temperature, or with the addition of 
excess Ir or Rh. Here we present evidence that this transition is driven by 
the lowering of the electronic density of states at the Fermi level of the alrV 
structure. The transition has long been thought to be second order, with 
a simple doubling of the LIq unit cell due to an unstable phonon at the R 
point (0||). We use first-principles calculations to show that all phonons at 
the R point are, in fact, stable, but do find a region of reciprocal space where 
the LIq structure has unstable (imaginary frequency) phonons. We use the 
frozen phonon method to examine two of these modes, relaxing the structures 
associated with the unstable phonon modes to obtain new structures which 
are lower in energy than LIq but still above alrV. We examine the phonon 
spectra of these structures as well, looking for instabilities, and find further 
instabilities, and more relaxed structures, all of which have energies above 
the cuIrV phase. In addition, we find that all of the relaxed structures, stable 
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and unstable, have a density comparable to the LIq phase (and less than 
the alrV phase), so that any transition from one of these structures to the 
ground state will have a volume change as well as an energy discontinuity. 
We conclude that the transition from LIq to alrV is probably weakly first 
order. 

Keywords: Structural Phase Transitions, Jahn- Teller, Electronic Structure, 
Density Functional Theory, Ordered Intermetallic Alloys, 



1. Introduction 

A major goal of computational condensed matter physics is the determi- 
nation of structural properties for compounds. Once the structure has been 
determined, other properties, e.g., strength, ductility, electronic properties 
(including superconductivity), magnetism, etc., can be determined. Several 
methods exist to attack this problem. One possibility is to exhaustively 
search an experimental database of known structures determining the low 
energy structures for each composition of the target materials. Other meth- 
ods use first-principles calculations on a small set of target structures to de- 
termine parameters which can be used to predict properties of more complex 
systems. These include the Cluster Expansion Method j2|, [s], 0, Isl, [oj, tight- 
binding parametrization methods 0, Isf, and atomistic potential methods j9[. 
Methods of the first type can be combined with methods of the second type 
for more through searches l3 • 



These programs, can, however, only be implemented with a thorough 
knowledge of the behavior of compounds. Exhaustive searches of a database 
require an extensive database to search, including all structures which are 
known, or thought likely, to form for a target system. Parameterized methods 
may not be able to reach all regions of phase space, and also need a database 



of structures for testing|lll. |12| . For both methods, then, it is useful to look 
at less-common structures found in nature. 

The alrV structure, shown in Figure [H is found in nature only in the pro- 
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totype compound |13l. |14| and its neighbor in the periodic table, RhV|15 
This orthorhombic structure, space group Cmmm (#65), can be viewed as 
a doubled unit cell distortion of the CsCl structure[l7| or the tetragonal LIq 
structure, space group P4/mmm (#123) [l^. Indeed, both the IrV and RhV 
phase diagrams |l9| show the LIq structure as the ground state for Vanadium- 
poor (a; = 40 — 48%) [Ir,Rh]i_2.Va:, and Vanadium-rich ctlrV is known to 
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Figure 1: The alrV structure. The box indicates the boundaries of the full orthorhombic 
unit cell. In the Llg or (CsCl) structure the Iridium atoms would be in a tetragonal (cubic) 
arrangement around each of the Vanadium atoms. 



transform to LIq at temperatures above 506° C 18 



The seemingly straightforward transition pathway from LIq to alrV was 
analyzed by Chen and Franzen 18| in the context of Landau theory: First, 
double the tetragonal LIq unit cell along the y and z directions, corresponding 
to a phonon at the R (0||) point in reciprocal space. Second, construct the 
primitive base-centered orthogonal unit cell, and allow the two Ir atoms in 
the cell to move opposite each other along the 2;-axis, while the V atoms 
move along the |/-axis. The resulting structure is symmetrically-equivalent 
to that of alrV and, assuming a second-order phase transition, relaxes into 
the alrV structure. 

This simple picture is not quite accurate. A true second order transi- 
tion, as described by Chen and Franzen, requires a continuous energy path 
from LIq to alrV. Thus the energy of the LIq structure would be lowered 
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by making an infinitesimal displacement of the type described above, cor- 
responding to a plionon instability at the R point. We have performed 
first-principles density functional theory calculations for the phonon frequen- 
cies at the R point, both by the frozen phonon method[2Q] and by linear 



response [21|, |22|, |23|. We find that all of the modes here have real frequen- 
cies, i.e., they do not lead to an instability. 

Thus the transition from the LIq phase to the alrV phase does not pro- 
ceed through the simple unit-cell doubling picture described above. Instead, 
as we shall show, the LIq phase has vibrational instabilities in another re- 
gion of the Brillouin zone. This paper will discuss the structures arising from 
these instabilities, and explore the possibility that a second-order transition 
might go through one of those phases. 

In Section [2] we show the results of density functional theory calculations 
for the energy and electronic structure of the alrV and LIq phases of IrV, 
RhV, and some neighboring compounds. We show that in IrV and RhV the 
LIq structure has a relatively high electronic density of states at the Fermi 



level, and so the transition to alrV has Jahn- Teller character [24 . 

In Section E] we look at the phonon spectra of the alrV and LIq phase. 
Not surprisingly, we find that the alrV phase has no phonon instabilities. 
We do, however, find that the LIq structure is vibrationally unstable on and 
near the line (x^^) in reciprocal space. 

In Section S] we use the frozen-phonon method on the unstable phonon 
modes to search for new structures which have lower energy than LIq. Such 
a search is computationally bound, so we looked at all unit cells with eight 
atoms or less, and two structures with 32 atom unit cells. In this range we find 
no instability which relaxes to the alrV phase, however it is not impossible 
that searching through larger unit cells would find such an instability, in 
which case the transition would indeed be second order. We do find several 
new structures, some of which have app arently never been seen in nature. 



and one which was previously known [16 



In Section [5] we discuss our results, and thoughts on the order of the 
Llo-alrV transition. 

2. Energetics and Electronic Structure of the LIq and alrV Phases 

All computations were made using the Kohn-Sham[25j| formulation of 



Density Functional Theory 26|| with the Perdew-Burke-Ernzerhof [27[] gener- 



alized gradient approximation. Depending on our needs, we used the Vienna 



4 




Figure 2: Energy /volume curves for IrV and RhV in the Llg and alrV structures, deter- 
mined from VASP and QE calculations as discussed in the text. For ease of comparison 
we have set the minimum energy of the LIq phase to zero. 



ab initio Simulation Package (VASP) [281, |29|, |30| with projector au gme nted 



wave (PAW) potentials [311], or the Quantum Espresso (QE) package |32[| with 
the supplied ultra-soft pseudopotentials for Ir, Rh, and V. 

We used rather large plane wave cutoffs of 350 eV in VASP and 540 eV in 
QE to ensure convergence. Similarly, we used 12 x 12 x 8 P-centered k-point 
meshes in both the LIq and alrV structures, leading to 140 and 215 points, 
respectively, in the irreducible Brillouin zones. We summed over electronic 
states using a Fermi-Dirac distribution 33(] with a temperature of 65 meV 
(0.005 Rydbergs). These values give well-converged results. 

We show our results for IrV and RhV in Table [1] and Fig. [2l In both 
cases we see that the VASP PAW potentials and the QE ultrasoft pseudopo- 
tentials are in excellent agreement with one another, and in good agreement 
with experiment, within the usual errors of Density Functional theory in the 
generalized gradient approximation. In both compounds the alrV state is 
below the LIq state by approximately 55 meV/formula unit, in agreement 
with experiment. 

From experiment 18j we know that LIq is the preferred high temperature 
structure for alrV. Using the COMSUBS routine from the ISOTROPYjsi, 
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package, we find that an orthorhombic distortion of the LIq structure 
reduces the symmetry from space group P4/mmm to Cmmm, the space 
group of the alrV structure. It is plausible to argue [l8| that this lower- 
ing of symmetry is the pathway for the LIq — )■ alrV phase transition. The 
symmetry-breaking character of the transition is evident from Fig. [3], where 
we plot the electronic density of states of both phases near the Fermi level. 
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Table 1: Equilibrium lattice constants (in A), atomic positions, and equilibrium bulk 
modulus (Kq, in GPa) for the LIq and alrV structures of IrV and RhV determined from 
experiment [13l [lit . VASP PAW calculations, and QE ultrasoft pseudopotential calcula- 
tions. Note that we give the primitive tetragonal lattice parameters for the LIq structure, 
rather than the more-common face-centered tetragonal setting, so that when c/a = 1 the 
LIq structure reduces to the cubic CsCl structure. The final row of the table shows the 
energy difference between the two phases (meV/formula unit). 



IrV 



RhV 



Llo Structure 





Exp. 


VASP 


QE 


Exp. 


VASP 


QE 


a 


2.749 


2.755 


2.762 


2.754 


2.739 


2.745 


c 


3.651 


3.677 


3.668 


3.599 


3.660 


3.648 


Ko 




272 


277 




230 




alrV Structure 


a 


5.791 


5.838 


5.816 


5.78 


5.849 


5.813 


b 


6.756 


6.759 


6.767 


6.65 


6.707 


6.725 


c 


2.796 


2.814 


2.823 


2.78 


2.792 


2.802 


Ko 




271 


276 


227 






Ir/Rh (4j) 


0.22 


0.216 


0.216 




0.214 


0.215 


V(4g) 


0.28 


0.296 


0.297 




0.296 


0.297 


A E 






53.6 


67.0 




50.4 


59.7 
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Figure 3: Electronic density of states for the LIq (red) and alrV (blue) phases of IrV 
(left) and RhV (right), found by smearing the eigenvalues computed by VASP using a 
Fermi-Dirac distribution at T = 65 meV (5 mRY) . ^33.] 



As we are only interested in the overall behavior of the density of states, we 
compute these curves by smearing out each eigenvalue found by VASP using 
a Fermi distribution at a temperature of 5 mRy. We see that at the Fermi 
level the density of states is twice as large in the LIq phase as it is in the 
alrV phase. This is consistent with a Jahn- Teller-like symmetry breaking 



and phase transition [2J|. As an aside, we also note that there is a minimum 
in the LIq density of states just above the Fermi level. This is consistent with 
the phase diagram of IrV(l9'], which shows that Iridium- rich IrV has the Llg 
structure. Assuming the additional Ir replaces V in the LIq structure, and 
using a rigid-band model, we can see how adding Ir to the system would raise 
the Fermi level and lower the density of states, leading to a more stable LIq 
phase. We find a similar, though less pronounced, state of affairs in RhV, as 
is also seen in Fig. El 

We investigated the possibility that the alrV structure might be found 
in compounds neighboring IrV and RhV in the periodic table. The results 
are shown in Table |2l Except for CoPt, all of the compounds shown here 
have 14 electrons in the valence band, just as in IrV and RhV. In all cases 
except CoPt we found that the alrV structure was locally stable, that is, 
if we started from the lattice parameters of the alrV, replacing Ir and V 
by the indicated atoms, and allowed the system to relax, it remained in the 
alrV structure and did not relax to the higher symmetry LIq structure. This 
suggests that the formation of the alrV structure is a result of a Fermi level 
effect associated with 14 electrons in a tetragonal unit cell. We did not check 
the elastic and vibrational stability of the resulting structure. However the 
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Table 2: Energy difference, in meV per formula unit between the LIq and alrV structures 
for the given compounds, as determined by VASP/PAW/PBE calculations. A positive 
number indicates that the Llo structure is lower in energy. An energy difference of zero 
indicates that the alrV structure relaxes into the LIq structure. Ground state structures 
are from the ASM Online Handbookfl^ . 

Compound CoPt CoV IrNb IrTa IrV RhNb RhTa RhV 



Ground State LIq D8b LIq 



alrV Llo 



alrV 



A E 







25.3 59.0 79.9 -53.6 50.1 70.6 -50.4 
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Figure 4: Electronic density of states of the ground-state LIq and hypothetical alrV phases 
of IrNb, found by smearing the eigenvalues computed by VASP. 

LIq structure was lower in energy for all compounds except IrV and RhV. 
Also, as shown in Fig.|H for these compounds the density of states of the LIq 
structure is very near to that of the alrV structure, precluding a Jahn- Teller 
lowering of the energy with decreased symmetry. We therefore conclude that 
it is unlikely to find the alrV structure anywhere except in IrV and RhV. 

3. Stability of the Llo and alrV Phases 

A system is stable, or metastable, only if it has real phonon frequen- 
cies for all k-points and modes except the acoustic phonons at the F point, 
which are guaranteed to vanish because of translational symmetry. Since the 
elastic constants are related to the long-wavelength behavior of the acoustic 
phonons[37] this implicitly includes the Born stability criteria [38|]. This is for- 
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Figure 5: The base-centered orthorhombic alrV structure. Left: Phonon frequencies IrV at 
equilibrium, determined using the hnear response method with ultrasoft pseudopotentials 
from the Quantum Espresso package. Right: First Brihouin zone. Symmetry hues are 
described by the notation of Miller and Love[36j. 



tunate, since the orthorhombic alrV structure would require us to compute 
nine elastic constants. 

To check the stability of the alrV phase and the possible meta-stability 
of the LIq phase we computed phonon frequencies throughout the respective 



Brillouin zones using linear response 21|, |2^, |23[ , as implemented in the Quan 



turn Espresso package PHon code[32|]. Phonon frequencies were computed on 
a reciprocal space grid ("q-points") using an 8 x 8 x 4 mesh for the LIq struc- 
ture (45 points in the irreducible Brillouin zone) and a 6 x 6 x 4 mesh for 
the alrV structure (39 points). Linear response calculations at these points 
yield a reciprocal-space dynamical matrix. These matrices are transformed 
into a real-space dynamical matrix which can be used to compute phonon 
frequencies over the entire Brillouin zone. 

This approach can lead to aliasing at points off of the q-point mesh. To 
check that this does not occur we also performed frozen-phonon calculations [2^ 
wherein we construct a supercell commensurate with the given q-point and 
measure the change in energy as a function of atomic displacement within the 
supercell. To do this we used the program FROZSL, part of the ISOTROPY 
package [s^. 35], with atomic displacements of 0.1 Bohr. Electronic structure 
calculations were then performed with the same energy cutoff as the original 
unit cell, and when possible the same k-point mesh, albeit folded back into 
the smaller Brillouin zone associated with the supercell. 

In the following, in interest of saving space, we only show the results for 
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Figure 6: The tetragonal LIq structure. Left: Phonon frequencies for IrV at its minimum 
energy structure, determined using the linear response method, using ultrasoft pseudopo- 
tentials from the Quantum Espresso package. Frequencies below zero (and in red) are 
unstable modes which actually have imaginary frequencies. The dots represent calcula- 
tions using the frozen phonon method with displacements computed by FROZSL. The 
second panel shows the frequencies along the line (a:;^^). Right: First Brillouin zone, with 
symmetry lines labeled according to the notation of Miller and Love[36i|. The arrow from 
U to S indicates the region of unstable phonons. 



IrV. Our calculations for RhV show a similar pattern. 

Fig. [5] shows the phonon frequencies of the equilibrium structure of alrV. 
High symmetry points and lines are labeled according to the convention of 
Miller and Love[36]. All the phonons have real frequencies, confirming that 
alrV is at least a metastable structure for IrV, as suspected. There is no 
evidence of any aliasing in these calculations. 

Figure [6] shows the phonon frequencies of the LIq phase of IrV, along 
symmetry lines labeled according to the convention of Miller and Love 36 . 
We show results of frozen phonon calculations, showing good agreement be- 
tween the two techniques, even off the linear response q-mesh. As expected, 
there are regions of reciprocal space with imaginary phonon frequencies. Sur- 
prisingly, this region is not near the R point (0||) as described by Chen and 



Franzen[18|, but along and near a line from the midpoint of the U line to the 
midpoint of the S line, as shown in more detail in the second panel of the 
figure. We will discuss the implications of this instability in the next section. 



4. Searching for Low Energy Structures 

The previous section showed that the alrV structure is the ground state 
(or, at least, a low-energy metastable state) of IrV, while the LIq structure is 
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unstable to phonons along and near the line (x^^) in the Brillouin zone, but 
not at the R point, where the phonon modes are real, invalidating the second- 
order phase transition scenario of Chen and Franzenfis!]. In this section we 
will look to see if it is possible to find another continuous transition path 
from LIq to alrV. 

In the frozen phonon method, we displace atoms from their equilibrium 
sites in a supercell consistent with the phonon wave vector and in directions 
which maintain the symmetry of the phonon mode. The frequency of the 
mode is then directly related to the square root of the curvature of the energy 
as a function of atom displacement. A mode with an imaginary frequency 
will then as a matter of course be related to a supercell calculation which 
has a negative curvature, leading to supercells with energy lower than the 
original state. These structures can then be relaxed, leading to new, or at 
least different, structures. 

As an example of this, consider the unstable LIq phonon at the point 
(iii) high symmetry "S" line. The linear response calculations find 

this phonon frequency to be 81t cm~^. To generate a supercell for this cal- 



culation we used the FROZSL code from the ISOTROPY package |3J, |35 
This package tells us that this point has phonons in three irreducible repre- 
sentations, each with two associated modes. For each representation we run 
three total energy calculations for displacements of the atoms in different 
directions, with the maximum displacement of 0.1 Bohr. The energy differ- 
ences between these structures and the ground state are used to determine 
the dynamical matrix at this point, and the mode frequencies. 

For the equilibrium LIq parameters found in Table [H the S3 irreducible 



representation (in the notation of Miller and Love[36|) has one mode with an 



imaginary frequency of 74 cm~^. (The discrepancy between linear response 
and frozen phonon calculations is due to the use of different k-point meshes 
and the anharmonicity of the mode. For our purposes we are only interested 
in showing that the mode is unstable in both cases, and so will not try to 
refine the calculations to improve the agreement between the two.) 

Diagonalizing the dynamical matrix allows us to find a supercell with 
displaced atoms which has an energy lower than the LIq phase. This supercell 
has space group Cmmm (#65), and with the appropriate choice of origin 
has Ir atoms at the (4e) and (4h) Wyckoff positions, and V atoms at the 
(2a), (2b), and (4j) positions. This is crystallographically equivalent to the 
intermetallic GaaPts structure jsoj, and we will refer to it as such. Upon 
relaxation we find a minimum energy structure which has approximately the 
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Figure 7: Linear-response phonon spectra of two eight-atom unit cells derived from un- 
stable phonon modes of LIq IrV. Left: The metastable RhsVa structure. Right: The 
Amm2 structure with Ir atoms at two (4d) Wyckoff positions and V atoms at four (2b) 
WyckofF positions. The imaginary frequencies of the unstable modes are shown as negative 
frequencies and highlighted in red. 



same density as, and an energy 20.5 meV/formula unit below, the relaxed 
LIq structure of Tabled] Note that this is still well above the energy of the 
alrV structure. 

In a similar fashion, if we look at the phonons at 0^|, on the "U" line, we 
find an imaginary U3 mode with a frequency of 70 cm~^. This mode leads to 
another supercell with space group Cmmm, but now the Ir atoms are on the 
(4e) and (4g) Wyckoff sites, while the V atoms occupy a pair of (4j) Wyckoff 
sites. The relaxed structure again has a density comparable to LIq, but its 
energy is only 11.0 meV/formula unit below the LIq structure minimum. In 
the discussion below will refer to this structure as Cmmm. 

Neither the GaaPts nor the Cmmm structure will relax to the ground 
state alrV structure. However, both structures have imaginary frequency 
long-wavelength optical mode phonons. If we venture away from the F point 
we will have to deal with frozen-phonon calculations on unit cells having 
hundreds of atoms, so we will only examine the unstable modes at F in both 
structures. 

First consider the imaginary frequency Fj mode of the Ga^Pt^ struc- 
ture. It is associated with a supercell of space group Amm2 (#38), with 
Ir atoms at the (4d) and (4e) Wyckoff sites and V atoms on pairs of (2a) 
and (2b) sites. This structure is crystallographically equivalent to the RhsVs 
structure (igI. 40 | and will be referred to by that name. The relaxed struc- 
ture is 41.1 meV/formula unit below the LIq structure, with similar density. 
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This structure is metastable, as seen in the hnear-response phonon spectrum 
plotted on the left-hand side of Fig. [71 

The second structure, derived from the unstable Fg mode of Cmmm, 
again with space group Amm2 (#38), puts the Ir atoms on two (4d) sites 
and the V atoms on four distinct (2b) sites. The relaxed structure has an 
energy 24.1 meV/formula unit below the Llo structure, again with a similar 
density. This structure has not been found in the intermetallic literature. It 
is vibrationally unstable, as can be seen from the phonon spectrum on the 
right-hand side of Fig. [71 

The four structures GasPts, RhsVs, Cmmm and Amm2 are the only 
structures with eight or fewer atoms in the primitive cell which can be gener- 
ated by relaxing the unstable phonon modes of the LIq structure. If we wish 
to examine larger unit cells, we can continue along this line of research ad in- 
finitum, or until we find a supercell where the atoms relax to the ground state 
alrV structure, a supercell with energy below the ground state structure, or 
we run out of computational resources. In most cases, including this one, we 
wiU reach the latter limit first. We did compute the relaxed energies of two 
supercells associated with unstable phonon modes in the Amm2 structure, 
and one with an unstable mode of the Cmmm structure. These structures 
had thirty-two atoms in the supercell, and energies below the parent struc- 
ture but above the metastable RhsVs structure. 

Tables [31 [Hand [SI summarize all of the calculations discussed here, includ- 
ing the structural derivation, space group, lattice constants, atomic positions, 
energy, and stability. 

5. Discussion 

We have shown by direct calculation that the transition from the LIq to 
the alrV structure of both IrV and RhV is a result of a Jahn- Teller driven 
distortion of the high-symmetry unit cell. This distortion does not result 
from a zone-doubling unstable phonon at the R point of the LIq Brillouin 
zone, and so the simple Landau theory picture does not hold. This does not 
completely eliminate the possibility that the transition is second order, for 
we found an entire region of reciprocal space, in and around the line (x^^), 
where the LIq structure has vibrational instabilities. We examined several 
of those instabilities, and found that they lead to numerous new structures, 
most of which have yet to be seen in ordered intermetallic systems and most 
of which are vibrationally unstable. The only phase which we found that 
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Table 3: Structural, energy, and vibrational stability results for the structures of IrV 
discussed in this paper. All calculations use the QE pseudopotentials discussed in the 
text. "Source" indicates the origin of the unit cell, either from experiment or an unstable 
phonon in the indicated structure. "Atoms" is the number of atoms in the primitive cell. 
Lattice constants are given in the standard crystallographic convention, e.g. the lattice 
constants of the full orthorhombic cell are given for the Cmmm and Cmcm space groups. All 
primitive cells have a = 7 = 90°. The rows "Ir" and "V" give the Wyckoff positions of 
the atoms. "Volume" is the minimum energy volume of the structure per formula unit, in 
(A^). "Energy" is the energy of the structure below the LIq structure, in meV/formula 
unit. 



Structure 


Llo 


Cmmm 


Cmcm 


Source 


Exp. 


U3(0i^) Llo 
phonon 


A3(00i) 

Cmmm phonon 


Space 
Group 


P4 / mmm 


Cmmm 


Cmcm 


Atoms 


2 


8 


32 


a (A) 


2.762 


7.260 


7.198 


b(A) 


2.762 


11.110 


11.166 


c(A) 


3.668 


2.774 


11.136 




90 


90 


90 


Ir 


(la) (000) 


(4e)(H0) 
(4g)(.267 0) 


(4d)(1^0) 
(4e)(.269 0) 
(4g)(.273 .254 \) 
(4g)(.268 .002 \) 


V 


(Id) (i|0) 


(4j)(0 .364 \) 
(4j)(0 .115 i) 


(41) (0 .358 .134) 
(4f)(0 .632 .117) 
(4f)(0 .123 .132) 
(4f)(0 .895 .117) 


Volume 


27.97 


27.97 


27.97 


Energy 





10.98 


20.61 


Stability 


Unstable 


Unstable 


Unknown 
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Table 4: Continuation of Table l3l 



Structure 




Amm2 


Imm2 


Source 


Ssdi^) Llo 
phonon 


Fg Cmmm 
phonon 


Hi(iOO) 
Amm2 phonon 


Space 
Group 


Cmmm 


Amm2 


Imm2 


Atoms 


8 


8 


32 


a (A) 


7.222 


2.788 


7.176 


b(A) 


7.800 


7.146 


11.144 


c(A) 


3.974 


11.230 


11.190 


/3 


90 


90 


90 


Ir 


(4e)(H0) 
(4h)(.273 i) 


(4d)(0 .230 .752) 
(4d)(0 .270 .998) 


(4c) (.256 .252) 
(8e)(.268 2.50 .252) 
(4c) (0.223 .752) 
(4c) (.256 .998) 
(8e)(.268 .250 .998) 
(4c)(.223 .498) 


V 


(2a) (000) 
(2b)(li0) 
(4j) (0 .218 i) 


(2b)(| .352) 
(2b)(i .625) 
(2b)(| .125) 
(2b)(i .899) 


(4d)(0 .131 .352) 
(4d)(0 .381 .355) 
(4d)(0 .132 .625) 
(4d)(0 .384 .625) 
(4d)(0 .130 .125) 
(4d)(0 .380 .125) 
(4d)(0 .131 .898) 
(4d)(0 .381 .895) 


Volume 


27.98 


27.97 


27.96 


Energy 


20.25 


24.15 


25.12 


Stability 


Unstable 


Unstable 


Unknown 
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Tabic 5: Continuation of Tables [5bnd 1^ 



Cm 
Ai(OiO) 
Amm2 
phonon 



Rh,V 



5 V3 



alrV 



GagPts 



phonon 



Exp. 



Cm 



Amm2 



Cmmm 



32 



13.303 
11.147 
6.652 
65.3 



4.070 
7.029 
7.822 

90 



5.816 
6.767 
2.823 
90 



J2a 
(2a 
(2a 
(2a 
(2a 
(2a 
(2a 
(2a 
(4b 
(4b 
(4b 
(4b 



.999 .504) 
.265 .974) 
.485 .534) 
.755 .994) 
.870 .256) 
.140 .715) 
.360 .276) 
.626 .745) 
.992 .249 .521) 
.260 .252 .983) 
.365 .748 .267) 
.133 .251 .729) 



"(4b 
(4b 
(4b 
(4b 
(4b 
(4b 
(4b 
(4b 



(4d)(0 .219 .724) 
(4e)(| .281 .001) 



(4j)(0 .216 i) 



.178 .130 .344) 
.175 .380 .362) 
.312 .134 .625) 
.312 .385 .625) 
.064 .125 .117) 
.060 .375 .133) 
.447 .130 .906) 
.450 .380 .888) 



(2a) (0 .032) 
(2a) (0 .457) 
(2b)(| .218) 
(2b)(| .793) 



(4g)(.297 0) 



27.97 



27.97 



27.78 



25.67 



41.14 



67.25 



Unknown 



Met ast able 



Ground State 
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is metastable is the experimentally observed RI15V3 phase, albeit in a 50-50 
composition. 

However, none of the unstable structures shows any sign of relaxing 
into the ground state structure. In fact, looking at Tables |3]|5], we see 
that all of the structures except alrV have approximately the same volume, 
27.97 A^/formula unit in IrV. The ground state alrV structure, one the other 
hand, has a somewhat smaller volume, 27.78 A^/formula unit. While this 
may seem insignificant, the volume change from one phase to another is a 
signal of a first-order transition, although only weakly first order here. 

This paper also shows a mechanism for generating new candidate in- 
termetallic phases: look for vibrationally unstable modes in a high energy 
structure, construct a supercell which will mimic that mode within the frozen- 
phonon calculation, and relax the cell. Using this method we found five new 
structures, as well as one, RhsVs, which had been seen before but which we 
had not been aware of until this research started. Many more new structures 
can undoubtedly be derived from just this system, but further research in this 
area is currently restricted by the time needed to search the Brillouin zone 
of a given crystal. For the eight atom supercells computing the phonon spec- 
trum at 45 q-points with a reasonable number (< 64) of processors took days 
to weeks to complete, depending on the symmetry of the crystal. Calculating 
the phonons for larger unit cells will of course take longer, although we will 
have fewer q-points to consider, alleviating some of the 0[A^^] increase due to 
the larger unit cell. One of our research goals will be to construct a "set-and- 
forget" mechanism which will search the Brillouin zone of an initial structure 
(here LIq) find all unstable modes with supercells containing a given number 
of atoms, relax those modes, and repeat, until all such structures have been 
found. Even restricting ourselves to all possible binary intermetallics[lj this 
will require a Grand Challenge computational program. 
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